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1. Preamble 



In the past two decades or so, we have learned how to understand crackling noise 
in a wide variety of systems. We review here the basic ideas and methods we 
use to understand crackling noise — critical phenomena, universality, the renor- 
malization group, power laws, and universal scaling functions. These methods 
and tools were originally developed to understand continuous phase transitions 
in thermal and disordered systems, and we also introduce these more traditional 
applications as illustrations of the basic ideas and phenomena. 

We focus largely on crackling noise in magnetic hysteresis, called Barkhausen 
noise. These lecture notes are distilled from a review article written with Karin 
Dahmen and Christopher Myers [1], from a book chapter written with Karin 
Dahmen and Olga Perkovic [2], and from a chapter in my textbook [3]. 

2. What is crackling noise? 

Many systems, when stressed or deformed slowly, respond with discrete events 
spanning a broad range of sizes. We call this crackling noise. The Earth crackles, 
as the tectonic plates rub past one another. The plates move in discrete earth- 
quakes (Fig.QJ), with many small earthquakes and only a few large ones. If the 
earthquake series is played as an audio clip, sped up by a factor of ten million, 
it sounds like crackling [4]. A piece of paper [6] or a candy wrapper [10] will 
crackle as it is crumpled (try it!), emitting sharp sound pulses as new creases 
form or the crease pattern reconfigures (Fig. 12a)). Paper also tears in a series 
of avalanches [7]; fracture in many other systems exhibit avalanche distributions 
and jerky motion [11, 12]. Foams (a head of beer) move in jerky avalanches as 
they are sheared [13], and as the bubbles pop [14]. Avalanches arise when fluids 
invade porous media in avalanches [15, 16] (such as water soaking into a sponge). 
I used to say that metals were an example of something that did not crackle when 
bent (permanently plastically deformed), but there is now excellent data that ice 
crackles when it is deformed [8], and recent data on micron-scale metal defor- 
mation also shows crackling noise [17, 18]. We will focus here on Barkhausen 
noise, the magnetic pulses emitted from (say) a piece of iron as it is placed inside 
an increasing external field (see Fig.^a)). 
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Energy radiated (10 joules) 




Fig. 1. Earthquake sizes, (a) Earthquake energy release versus time in 1995. There are only a 
few large earthquakes, and many small ones. This time series, when sped up, sounds like crackling 
noise [4]. (b) Histogram of the number of earthquakes in 1995 as a function of their magnitude M . 
Notice the logarithmic scales; the smallest earthquakes shown are a million times smaller and a thou- 
sand times more probable than the largest earthquakes. The fact that the earthquake size distribution 
is well described by a power is the Gutenberg-Richter law [5]. 



All of these systems share certain common features. They all have many more 
tiny events than large ones, typically with a power-law probability distribution 
(Figs Ufa), Ha)). For example, the histogram of the number of earthquakes of a 
given size yields a straight line on a log-log plot (FigureQJ). This implies that the 
probability of a large earthquake goes as its size (energy radiated) to a power. The 
tearing avalanches in paper, dislocation avalanches in deformed ice, and rupture 
avalanches in steel all have power-law avalanche size distributions(Fig. 12b)). 
Because crackling noise exhibits simple, emergent features (like these power law 
size distributions), it encourages us to expect that a theoretical understanding 
might be possible. 

What features of, say, an earthquake fault do we expect to be important for 
such a theoretical model? If earthquake faults slipped only in large, snapping 
events we would expect to need to know the shape of the tectonic plates in order 
to describe them. If they slid more-or-less smoothly we anticipate that the nature 
of the internal rubble and dirt (fault gouge) would be important. But since earth- 
quakes come in all sizes, we expect that neither the microscopic rock-scale nor 
the macroscopic continental-scale details can be crucial. What, then, is important 
to get right in a model? 

An important hint is provided by looking at the dynamics of an individual 
avalanche. In many (but not all) systems exhibiting crackling noise, the avalanches 
themselves have complex, internal structures. Fig. |4] shows that the avalanches 
producing individual Barkhausen pulses in magnets proceed in an irregular, jerky 
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Fig. 2. (a) Crackling noise in paper V(t) for multiple pulses, from [6]. Again, there are many small 
pulses, and only a few large events, (b) Power laws of event sizes for paper tearing [7], plastic flow 
in ice [8], and brittle and ductile fracture in steel [9]. 



fashion, almost stopping several times in between sub-avalanches. Fig. |3Ja) 
shows that the spatial structure of an avalanche in a model magnet is also ir- 
regular. Both the time and spatial structures are self-similar; if we take one of 
the sub-avalanches and blow it up, it looks statistically much like the original 
avalanche. Large avalanches are built up from multiple, similar pieces, with each 
sub-avalanche triggering the next. It is the way in which one sub-avalanche trig- 
gers the next that is crucial for a theoretical model to get right. 

Our focus in these lecture notes will be on understanding the emergent be- 
havior in crackling noise (power laws and scaling functions) by exploring the 
consequences of this self-similar structure. In section [3] we illustrate the pro- 
cess by which models are designed and tested by experiments using Barkhausen 
noise. In section|4]we introduce the renormalization group, and use it to explain 
universality and self-similarity in these (and other) systems. In section|5]we use 
the renormalization group to explain the power laws characteristic of crackling 
noise, and also use it to derive the far more powerful universal scaling functions . 



3. Hysteresis and Barkhausen noise in magnets 

Microscopically, iron at room temperature is always magnetized; non-magnetic 
bulk iron is composed of tiny magnetic domains, whose north poles point in 
different directions giving a net zero magnetization. An external magnetic field, 
as in Fig. 0a), will attract a piece of iron by temporarily moving the domain 
walls to enlarge the regions whose north poles are aligned towards the south pole 
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Fig. 3. (a) Barkhausen noise experiment. By increasing an external magnetic field H(i) (bar 
magnet approaching), the magnetic domains in a slab of iron flip over to align with the external field. 
The resulting magnetic field jumps can be turned into an electrical signal with an inductive coil, and 
then listened to with an ordinary loudspeaker. Barkhausen noise from our computer experiments can 
be heard on the Internet [4]. (b) Cross section of all avalanches in a billion-domain simulation of 
our model for Barkhausen noise at the critical disorder [1]. The white background is the infinite, 
spanning avalanche. 



of the external field. When the external field is removed, these domain walls do 
not completely return to their original positions, and the iron will end up partially 
magnetized; this history dependence is called hysteresis. The dependence of the 
magnetization of the iron M on the history of the external field H (the hysteresis 
loop) can have an interesting hierarchical structure of subloops (Fig. [6]). 

These hysteresis loops may look smooth, but when Fig.|6]is examined in detail 
we see that the magnetization grows in discrete jumps, or avalanches. These 
avalanches are the origin of magnetic Barkhausen noise. In section [XT1 we will 
develop a simple model for this Barkhausen noise. In section [J!2l we shall observe 
that our model, while capturing some of the right behavior, is not the correct 
model for real magnets, and introduce briefly more realistic models that do appear 
to capture rather precisely the correct behaviorQ 

3.1. Dynamical Random-field Ising model 

We introduce here a caricature of a magnetic material. We model the iron as a 
cubic grid of magnetic domains Si, whose north pole is either pointing upward 
(Si — +1) or downward (Si — —1). The external field pushes on our domain 
with a force H(t), which starts pointing down (H (t = 0) 0) and will increase 



'That is, we believe they are in the right universality class, section pTTl 
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Fig. 4. Internal avalanche time series (a) during a single experimental Barkhausen noise pulse [19], 
and (b) during a single avalanche in a simulation of magnetic Barkhausen noise [1]. The experiment 
measures the voltage V(t), which in this experiment measures the volume swept out per unit time 
by the moving magnetic domain wall; the simulation measures the number of domains flipped per 
unit time. In these two cases, the total area under the curve gives the size S of the avalanche. Notice 
how the avalanches almost stop several times; if the forcing were slightly smaller, the large avalanche 
would have broken up into two or three smaller ones. The fact that the forcing is just large enough to 
on average keep the avalanche growing is the cause of the self-similarity; a partial avalanche of size 
S will on average trigger one other of size S. 



with time. Iron is magnetic because a domain has lower energy when it is mag- 
netized in the same direction as its neighbors; the force on site i from the six 
neighboring domains Sj in our model is of strength J Y^j Sj . Finally, we model 
the effects of impurities, randomness in the domain shapes, and other kinds of 
disorder by introducing a random fielc0 hi, different for each domain and chosen 
from a normal distribution with standard deviation Rx J. The net force on Si is 
thus 

Ft = H (t) + J Sj + K (3.1) 

3 

corresponding to the energy function or Hamiltonian 

H = -jJ^SiSj - 5^(JT(t) + h)Si. (3.2) 

(ij) * 



2 Most disorder in magnets is not microscopically well described by random fields, but better 
modeled using random anisotropy or random bonds which do not break time-reversal invariance. For 
models of hysteresis, time-reversal symmetry is already broken by the external field, and all three 
types of randomness are probably in the same universality class. That is, random field models will 
have the same statistical behavior as more realistic models, at least for large avalanches and long 
times (see section |4~TV 
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Fig. 5. (a) Fractal spatial structure of an avalanche. This moderate-sized avalanche contains 
282,785 domains (or spins) [1]. The shading depicts the time evolution: the avalanche started in the 
dark region in the back, and the last domains to flip are in the upper, front region. The sharp changes 
in shading are real, and represent sub-avalanches separated by times where the avalanche almost 
stops (see Fig.|4). (b) Avalanche propagation. The avalanches in many models of crackling noise 
are first nucleated when the external stress induces a single site to transform. In this two-dimensional 
model, site #13 is the nucleating site which is first pushed over. The coupling between site #13 and 
its neighbors triggers some of them to flip (#4, 8, and 12), which in turn trigger another shell of 
neighbors (#15, 19, 7, 11, and 17), ending eventually in a final shell which happens not to trigger 
further sites (#6, 10, 18, and 20). The time series V(t) plotted in figure|4lb) is just the number of 
domains flipped in shell #t for that avalanche. 



This model is called the random-field Ising model (RFIM), and its thermal equi- 
librium properties have historically been studied as an archetypal disordered sys- 
tem with glassy dynamics. Here we ignore temperature (specializing to magnets 
where the thermal fluctuations are too small to de-pin the domain walls), and as- 
sume that each domain Si reorients whenever its local field F t changes sign, to 
minimize its energy. 

In our model, there are two situations where a domain will flip. It may be 
induced to flip directly by a change in the external field H(t), or it may be trig- 
gered to flip by the flipping of one of its neighboring domains. The first case 
corresponds to nucleating a new avalanche; the second propagates the avalanche 
outward from the triggering domain, see Fig. [2b). We assume that the external 
field H (i) changes slowly enough that each avalanche finishes (even ones which 
sweep over the entire magnet) before the field changes appreciably. 

These avalanches can become enormous, as in Fig.Qa). How can an avalanche 
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Fig. 6. Hysteresis loop and subloops: The magnetization in our model [20], as the external field H 
is ramped up and down. Our focus will primarily be on the upper, outer loop as the external field is 
ramped from — oo to oo. Barkhausen jumps (exploded region): The hysteresis loop appears smooth, 
but when examined in detail is composed of discrete, abrupt jumps in magnetization, corresponding 
to avalanches in the positions of the walls of the magnetic domains. This jerky magnetization is what 
emits Barkhausen noise. 



grow to over 10 5 domains, but then halt? In our model, it happens only when the 
disorder R and the field H(t) are near a critical point (Fig. [7). For large dis- 
order compared to the coupling between domains R 3> J = 1, each domain 
turns over almost independently (roughly when the external field H(t) cancels 
the local random field hi, with the alignment of the neighbors shifting the transi- 
tion only slightly); all avalanches will be tiny. For small disorder compared to the 
coupling R <C J, there will typically be one large avalanche that flips most of the 
domains in the sample. (The external field H (t) needed to nucleate the first do- 
main flip will be large, to counteract the field —6 J from the unflipped neighbors 
without substantial assistance from the random field; at this large external field 
most of the neighboring domains will be triggered, . . . eventually flipping most 
of the domains in the entire system). There is a critical disorder R c separating 
two qualitative regimes of behavior — a phase R > R c where all avalanches are 
small and a phase R < R c where one 'infinite' avalanche flips a finite fraction of 
all domains in the system (Fig. [7j. Neai0 R c , a growing avalanche doesn't know 
whether it should grow forever or halt while small, so a distribution of avalanches 
of all sizes might be expected — giving us crackling noise. 



3 One must also be near the field H (t) = H c where the infinite avalanche line ends. 
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Fig. 7. Phase transition in hysteresis and avalanche model; (a) one enormous avalanche (background) 
for small disorder R < R c , (b) avalanches at all scales at R = R c (c) many small avalanches for 
large disorder R > R c . (d) The upper branch of the hysteresis loop develops a macroscopic jump in 
magnetization below R c , whose size diverges as AAf oc (R — R c )@; at R = R c the magnetization 
has a power-law singularity M — M c ~ (H — Hc) 1 ^ . 

3.2. Real Barkhausen noise 

Is our model a correct description of Barkhausen noise in real magnets? Our 
model does capture the qualitative physics rather well; near R c it exhibits crack- 
ling noise with a broad distribution of avalanche sizes (Fig. |5Jb)), and the indi- 
vidual avalanches have the same kind of internal irregular dynamical structure as 
seen in real Barkhausen avalanches (Fig. |4j. 

However, when examined in detail we find that there are problems with the 
model. A typical problem might be that the theoretically predicted power-law is 
wrong (Fig. [8). The probability density V(S) of an avalanche of size S experi- 
mentally [21] decays as a power law D(S) ~ S~ T , with values of the exponent 
that cluster around either r = 1.27 or r = 1.5. In our model at R c and near 
H c we find r = 1.6 ± 0.06, which is marginally compatible with the data for 
the second group of samples (see Fig. [151' a)). If the power law were significantly 
different, our model would have been ruled out. We shall see in section[5]that the 
power laws are key universal predictions of the theory, and if they do not agree 
with the experiment then the theory is missing something crucial. 

However, there are two other closely related avalanche models that describe 
the two groups of experiments well. One is a front propagation or fluid invasion 
model (which preceded ours [15, 22-24]), also usually written as a random-field 
Ising model but here starting with an existing front and allowing new avalanches 
only contiguous to previously flipped domains (water can invade new pores only 
next to currently wet pores). The front-propagation model does well in describ- 
ing the dynamics and size distribution of the materials with r ~ 1.27. The other 
model can either be viewed as the motion of a flat domain wall [25, 26] or the 
mean-field limit of the random-field Ising model [27,28]; it has a value r = 3/2 
nicely describing the other class of materials. These models incorporate the long- 




Fig. 8. Avalanche size distributions for magnetic Barkhausen noise in various materials, from [21]. 
The materials fall into two families, with power law exponent D(S) ~ S~ T with r ~ 1.5 and 
t ~ 1.27. The cutoffs in the avalanche size distribution at large S are due to demagnetization 
effects, and scale as 5 m ax 

fe-V^fc (sectionEm Fig.[T7l 



range magnetic interactions between flipped domains, which experimentally pro- 
duce demagnetizing fields that induce the magnetic domain wall to move rather 
rigidly. 
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Fig. 9. Domain structure in a real magnet, measured with a scanning electron microscope, digitally 
processed to detect edges. The horizontal lines are magnetic domain walls. Courtesy of K. Zaveta, 
from [29]. 

It is easy to measure the exponents of a power-law decay, but it is difficult to 
measure them well. We could argue that r = 1.6 is close enough to some of 
the experimental measurements that our model is not ruled out. But even a brief 
investigation of the qualitative behavior in these systems leads us to abandon 
hope. Fig.|9]shows the domain-wall structure in a real magnet. The domain walls 
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are indeed flat, and the way in which they advance (looking at animated versions 
of this figure, or observing directly through a microscope) is not by nucleating 
new walls, but by motion of existing boundaries. 



4. Why Crackling Noise? 

In this section we introduce universality, the renormalization group, and self- 
similarity. These three notions grew out of the theory of continuous phase transi- 
tions in thermal statistical mechanics, and are the basis of our understanding of a 
broad variety of systems, from quantum transitions (insulator to superconductor) 
to dynamical systems (the onset of chaos) [3, chapter 12]. We will illustrate and 
explain these three topics not only with our model for crackling noise in mag- 
nets, but also using classic problems in statistical mechanics — percolation and 
the liquid-gas transition. The renormalization group is our tool for understanding 
crackling noise. 

4.1. Universality 

Consider holding a sheet of paper by one edge, and sequentially punching holes 
of a fixed size at random positions. If we punch out only a few holes, the paper 
will remain intact, but if the density of holes approaches one the paper will fall 
apart into small shreds. There is a phase transition somewhere in between where 
the paper first ceases to hold together. This transition is called percolation. 

Fig. [10] shows two somewhat different microscopic realizations of this prob- 
lem. On the top we see a square mesh of bonds, where we remove all but a 
fraction p of the bonds. On the bottom we see a lattice of hexagonal regions 
punched out at random. Microscopically, these two processes seem completely 
different (left figures). But if we observe large simulations and examine them at 
long length scales, the two microscopic realizations yield statistically identical 
types of percolation clusters (right FigsfTOl. 

Universality arises when the statistical morphology of a system at a phase tran- 
sition is largely independent of the microscopic details of the system, depending 
only on the type (or universality class) of the transition. This should not come as 
a complete surprise; we describe most liquids with the same continuum laws (the 
Navier-Stokes equations) despite their different molecular makeups, with only 
the mass density and the viscosity depending on microscopic details. Similarly 
most solids obey elasticity theory on scales large compared to their constituent 
particles. For these phases we understand the material independence of the con- 
stituent equations by taking a continuum limit, assuming all behavior is slowly 
varying on the scale of the particles. At critical points, because our system is 
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Fig. 10. Universality in percolation. Universality suggests that the entire morphology of the percola- 
tion cluster at p c should be independent of microscopic details. On the top, we have bond percolation, 
where the bonds connecting nodes on a square lattice are occupied at random with probability p; the 
top right shows the infinite cluster on a 1024 X 1024 lattice at the critical point. On the bottom, we 
have site percolation on a triangular lattice, where it is the hexagonal sites that are occupied. Even 
though the microscopic lattices and occupation rules are completely different, the resulting clusters 
look statistically identical at their critical points. (One should note that the site percolation cluster is 
slightly less dark. Universality holds up to overall scale changes, here up to a change in the density.) 
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rugged and complicated all the way down to the lattice scale, we will need a 
more sophisticated way of understanding universality ( section l4~2l >. 

Most thermal phase transitions are not continuous, and do not have structure 
on all scales. Water is water until it is heated past the boiling point, after which it 
abruptly turns to vapor (with a big drop in density p). Even boiling water, though, 
can show big fluctuations at high enough pressure. As we raise the pressure the 
boiling temperature gets larger and the density drop gets smaller, until at a certain 
pressure and temperature (P c , T c ) the two densities become equal (to p c ). At this 
critical point the H2O molecules do not know which phase they want to be in, 
and they show fluctuations on all scales. 

Figure [TTT a) shows that the approach to this critical point has universal fea- 
tures. If we divide the density by p c and the temperature by T c , many different 
liquids share the same p, T phase diagram, 



This is partly for mundane reasons; most of the molecules are roughly spherical, 
and have weak interactions. 

But figure fTTT b) shows the phase diagram of a completely different system, a 
magnet whose magnetization M(T) vanishes as the temperature is raised past its 
T c , The curve here does not agree with the curves in the liquid-gas collapse, but 
it can agree if we shear the axes slightly: 



Nature has decided that the 'true' natural vertical coordinate for the liquid-gas 
transition is not T but the line T = (p— A2) /A3. Careful experiments, measuring 
a wide variety of quantities, show that the liquid-gas transition and this type of 
magnet share many properties, except for smooth changes of coordinates like 
that in eq !4.2l Indeed, they share these properties also with a theoretical model — 
the (thermal, non-random, three-dimensional) Ising model. Universality not only 
tied disparate experiments together, it also allows our theories to work. 

The agreement between Figs [TTT a) and (b) may not seem so exciting; both 
are pretty smooth curves. But notice that they do not have parabolic tops (as 
one would expect from the maxima of a typical function). Careful measurements 
show that the magnet, the liquid-gas critical point, and the Ising model all vary as 
(1 — T/T c )° with the same, probably irrational exponent (3; the best theoretical 
estimates have (3 — 0.325 ± 0.005 [32, chapter 28]. This characteristic power 
law represents the effects of the large fluctuations at the critical points (the peaks 
of these graphs). Like the avalanche size distribution exponent r, /3 is a universal 
critical exponent (section[3T|l. 



p (T) = Ap (BT). 



(4.1) 



p (T) = A X M{BT) + A 2 + A 3 T. 



(4.2) 
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Fig. 1 1 . Universality, (a) Universality at the liquid-gas critical point. The liquid-gas coexistence 
lines (p(T)/ p c versus T/T c ) for a variety of atoms and small molecules, near their critical points 
(T c , p c ) [30]. The gas phases lie on the upper left; the liquid phase region is to the upper right; 
densities in between will separate into a portion each of coexisting liquid and gas. The curve is a 
fit to the argon data, p/p c = 1 + s(l - T/T c ) ± p (l - T/T C )P with s = 0.75, p = 1.75, 
and = 1/3 [30]. (b) Universality: ferromagnetic-paramagnetic critical point. Magnetization 
versus temperature for a uniaxial antiferromagnet MnF2 [31]. The tilted line in (a) corresponds to the 
vertical axis in (b). 
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4.2. Renormalization group 

Our explanation for universality, and our theoretical framework for studying 
crackling noise, is the renormalization group. The renormalization group starts 
with a remarkable abstraction: it works in an enormous 'system space'. Different 
points in system space represent different materials under different experimental 
conditions, or different theoretical models with different interactions and evolu- 
tion rules. For example, in Fig. [T2l' a) we consider the space of all possible mod- 
els and experiments on hysteresis and avalanches, with a different dimension for 
each possible parameter (disorder R, coupling J, next-neighbor coupling, . . . ) 
and for each parameter in an experiment (chemical composition, annealing time, 
. . . ). Our theoretical model will traverse a line in this infinite-dimensional space 
as the disorder R is varied. 

The renormalization group studies the way in which system space maps into 
itself under coarse-graining. The coarse-graining operation shrinks the system 
and removes microscopic degrees of freedom. Ignoring the microscopic de- 
grees of freedom yields a new physical system with the same properties at long 
length scales, but with different (renormalized) values of the parameters. As 
an example, figure Qj] shows a real-space renormalization-group 'majority rule' 
coarse-graining procedure applied to the Ising model. Many approximate meth- 
ods have been devised to implement this coarse-graining operation (real-space, 
e-expansions, Monte Carlo. . . ) which we will not discuss here. 

Under coarse-graining, we often find a fixed point S* in system space. All 
of the systems that flow into this fixed point under coarse-graining will share 
the same long-wavelength properties, and will hence be in the same universality 
class. 

Fig. [12 a) shows the case of a fixed point S* with one unstable direction. 
Points deviating from S* in this direction will flow away from it under coarse- 
graining. There is a surface C of points which do flow into the fixed point, 
which separates system space into two different phases (say, one with all small 
avalanches and one with one system-spanning, 'infinite' avalanche). The set C 
represents a universality class of systems at their critical points. Thus, fixed 
points with one unstable direction represent phase transitions. 

Cases like the liquid-gas transition with two tuning parameters (T c , p c ) deter- 
mining the critical point will have two unstable directions. What happens when 
we have no unstable directions? The fixed-point S a in Fig. [T^b) represents an 
entire region in system space sharing the same long-wavelength properties; it 
represents a phase of the system. Usually phases do not show structure on all 
scales, but some cases (like random walks) show this generic scale invariance. 

Sometimes the external conditions acting on a system naturally drive it to stay 
at or near a critical point, allowing one to spontaneously observe fluctuations on 
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Fig. 12. Renormalization-group flows [1]. (a) Flows describing a phase transition. The 

renormalization-group uses coarse-graining to longer length scales to produce a mapping from the 
space of physical systems into itself. Consider the space of all possible systems exhibiting mag- 
netic hysteresis (including, both real models and experimental systems). Each model can be coarse- 
grained, removing some fraction of the microscopic domains and introducing more complex dynam- 
ical rules so that the remaining domains still flip over at the same external fields. This defines a 
mapping of our space of models into itself. A fixed point S* in this space will be self-similar: be- 
cause it maps into itself upon coarse-graining, it must have the same behavior on different length 
scales. Points that flow into S* under coarse-graining share this self-similar behavior on sufficiently 
long length scales: they all share the same universality class, (b) Attracting fixed point [1]. Often 
there will be fixed points that attract in all directions. These fixed points describe typical behavior: 
phases rather than phase transitions. Most phases are rather boring on long length scales. In more in- 
teresting cases, like random walks, systems can exhibit self-similarity and power laws without special 
tuning of parameters. This is called generic scale invariance. (c,d) Flows for a front-propagation 
model [1]. The front propagation model has a critical field F c at which the front changes from a 
pinned to sliding state, (c) Coarse-graining defines a flow on the space of earthquake models. The 
critical manifold C, consisting of models which flow into S* * „, separating stuck faults from faults 
which slide forward with an average velocity v(F). (d) The velocity varies with the external force 
F across the fault as a power law v(F) ~ (F — F c )@. Clever experiments, or long-range fields, 
can act to control not the external field, but the position: changing the front displacement slowly sets 
v fa 0, thus self-tuning F fa F c . This is one example of self-organized criticality. 
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Fig. 13. Ising model at T c : coarse-graining. Coarse-graining of a snapshot of the two-dimensional 
Ising model at its critical point. Each coarse-graining operation changes the length scale by a factor 
B = 3. Each coarse-grained spin points in the direction given by the majority of the nine fine-grained 
spins it replaces. This type of coarse-graining is the basic operation of the real-space renormalization 
group. 
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all scales. A good example is provided by certain models of earthquake fault 
dynamics. Fig. [T2l c) schematically shows the renormalization-group flow for an 
earthquake model. It has a depinning fixed point representing the external force 
F across the fault at which the fault begins to slip; for F > F c the fault will 
slide with a velocity v(F) ~ (F — Fc) 13 . Only near F c will the system exhibit 
earthquakes of all sizes. Continental plates, however, do not impose constant 
forces on the faults between them; they impose a rather small constant velocity 
(on the order of centimeters per year, much slower than the typical fault speed 
during an earthquake). As illustrated in Fig. Q~2d), this naturally sets the fault 
almost precisely at its critical point, tuning the system to its phase transition. 
This is called self-organized criticality. 

In broad, the existence of fixed points under renormalization group coarse- 
graining is our fundamental explanation for universality. Real magnets, Ising 
models, and liquid-gas systems all apparently flow to the same fixed points under 
coarse graining, and hence show the same characteristic fluctuations and behav- 
iors at long length scales. 



5. Self-similarity and its consequences 

The most striking feature of crackling noise and other critical systems is self- 
similarity, or scale invariance. We can see this vividly in the patterns observed 
in the Ising model (upper left, Fig.[T3l), percolation (Fig.fTOb. and our RFIM for 
crackling noise (Fig. [5^). Each shows roughness, irregularities, and holes on all 
scales at the critical point. This roughness and fractal-looking structure stems 
at root from a hidden symmetry in the problem: these systems are (statistically) 
invariant under a change in length scale. 

Consider Fig. [HI depicting the self-similarity of the avalanches in our RFIM 
simulation at the critical disorder. The upper-right figure shows the entire system, 
and each succeeding picture zooms in by another factor of two. Each zoomed- 
in picture has a black 'background' showing the largest avalanche spanning the 
entire system, and a variety of smaller avalanches of various sizes. If you blur 
your eyes a bit, the figures should look roughly alike. This rescaling and blurring 
process is the renormalization-group coarse-graining transformation. 

How does the renormalization group explain self-similarity? The fixed point 
S* under the renormalization group is the same after coarse-graining (that's what 
it means to be a fixed point). Any other system that flows to S* under coarse 
graining will also look self-similar (except on the microscopic scales that are 
removed in the first few steps of coarse-graining, during the flow to S*). Hence 
systems at their critical points naturally exhibit self-similarity. 
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Fig. 14. Avalanches: scale invariance. Magnifications of a cross-section of all the avalanches in a 
run of our hysteresis model each one the lower right-hand quarter of the previous. The system started 
with a billion domains (1000 3 ). Each avalanche is shown in a different shade. Again, the larger scales 
look statistically the same. 
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This scale invariance can be thought of as an emergent symmetry, invariance 
under changes of length scale. In a system invariant under a translational sym- 
metry, the expectation of any function of two positions x\ and X2 can be written 
in terms of the separation between the two points, (g(xi,X2)) = G{%2 — Xi), 
In just the same way, scale invariance will allow us to write functions of N vari- 
ables in terms of scaling functions of N — 1 variables — except that these scaling 
functions are typically multiplied by power laws in one of these variables. 

5.1. Power laws 

Let us begin with the case of functions of one variable. Consider the avalanche 
size distribution D(S) for a model, say the real earthquakes in Fig.[Tfa), or our 
model for hysteresis at its critical point. Imagine taking the same system, but 
increasing the units of length with which we measure the system — stepping back, 
blurring our eyes, and looking at the system on a coarse-grained level. Let us 
multiply the spacing between markings on our rulers by a small amount B = 
(1 + e). After coarsening, any length scales in the problem (like the spatial extent 
L of a particular avalanche) will be divided by B. The avalanche size (volume) S 
after coarse-graining will also be smaller by some factoiQc = (1 + ce). Finally, 
the overall scale of D(S) after coarse-graining will be changed by some factor 
A = 1 + ae0 Hence under the coarse-graining we have 

L' = LJB = L/(l + e), 

S' = S/C = 5/(1 + ce), (5.1) 
D' = AD = D(l + ae). 

Now the probability that the coarse-grained system has an avalanche of size S' 
is given by the rescaled probability that the original system had an avalanche of 
size S = (1 + ce)S': 

D'(S r ) = AD{CS') = (1 + ae)D ((1 + ce)S')) . (5.2) 

Here D'(S') is the distribution measured with the new ruler: a smaller avalanche 
with a larger probability density. Because we are at a self-similar critical point, 
the coarse-grained distribution D'(S') should equal D(S'). Making e infinitesi- 

4 If the size of the avalanche were the cube of its length, then c would equal three since (1+e) 3 = 
1 + 3e + 0(t 2 ). In general, c is the fractal dimension of the avalanche. 

5 The same avalanches occur independent of your measuring instruments! But the probability 
D(S) changes, because the fraction of large avalanches depends upon how many small avalanches 
you measure, and because the fraction per unit S changes as the scale of S changes. 
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mal leads us to a differential equation: 



D(S>) 







D'(S') = (l + ae)D((l + ce)S'), 
aeD + ceS'-t-;, 



dD 

lis 



aD 



(5.3) 



which has the general solution 



D = D S- a/c . 



(5.4) 



Because the properties shared in a universality class only hold up to over- 
all scales, the constant Dq is system dependent. However, the exponents a, c, 
and a/c are universal — independent of experiment (within the universality class). 
Some of these exponents have standard names: the exponent c giving the frac- 
tal dimension of the avalanche is usually called df or \jav. The exponent a/c 
giving the size distribution law is called r in percolation and in most models 
of avalanches in magnets^] and is related to the Gutenberg-Richter exponent for 
earthquakes^ (Fig. |TJb)). Most measured quantities depending on one variable 
will have similar power-law singularities at the critical point. For example, the 
distribution of avalanche durations and peak heights also have power-law forms. 
This is because power laws are the only self-similar functions. If f(x) = x~ a , 
then on a new scale multiplying x by B, f(Bx) = B~ a x~ a oc f(x). 

Crackling noise involves several power laws. We've seen that the probability 
of having an avalanche of size S goes as S~ T . In our RFIM model for hysteresis, 
if one is at a distance R — R c from the critical point, there will be a cutoff in the 
avalanche size distribution. The typical largest spatial extent L of an avalanche is 
called the correlation length £ , which scales as £ ~ (R—R c )~ v . The cutoff in the 
avalanche size S scales as (R—R c )~ a (FigfTTTi. In other models, demagnetization 
effects (parameterized by k) lead to cutoffs in the avalanche size distribution [21], 
with analogous critical exponents v% and o^. The size and spatial extent of the 
avalanches are related to one another by a power law S ~ L d f , where dt = 1 /<rv 
is called the fractal dimension. The duration of an avalanche goes as L z . The 
probability of having an avalanche of duration T goes as T~ a , where a = (r — 



6 In our RFIM for hysteresis, we use r to denote the avalanche size law at the critical field and 
disorder (D(S, R c , H c ) ~ S~ T ); integrated over the hysteresis loop D; nt oc S~ T with f = 
r + a/38. 

7 We must not pretend that we have found the final explanation for the Gutenberg-Richter law. 
There are many different models that give exponents 2/3, but it remains controversial which of 
these, if any, are correct for real-world earthquakes. 
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Fig. 15. (a) Universal Critical Exponents in Various Spatial Dimensions. [1] We test our e- 
expansion predictions [34] by measuring the various critical exponents numerically in up to five 
spatial dimensions [35, 36]. The various exponents are described in the text. All of the exponents are 
calculated only to linear order in e, except for the correlation length exponent u, where we use results 
from other models. The agreement even in three dimensions is remarkably good, considering that 
we're expanding in 6 — D where D = 3! (b) Universal Critical Exponents vs. Experiment [21]. 
The exponent r gives the probability density ~ S~ T of having an avalanche of size S, a gives the 
probability density ~ T~ a of having an avalanche of duration T, and the exponents cry. and 
describe the cutoff in the avalanche sizes and durations, respectively, due to the demagnetizing field 
(see [21]). The experimental samples group naturally into two groups, described by the mean-field 
and front-propagation universality classes; none appear to be described well by our nucleated RFIM. 



Vj/avz + 1@ In our RFIM, the jump in magnetization goes as (i? — i? c ) /3 , and 
at R c the magnetization (M - M c ) - (H - H c f/ S (Fig.^d)). 

To specialists in critical phenomena, these exponents are central; whole con- 
versations often rotate around various combinations of Greek letters. We know 
how to calculate critical exponents from various analytical approaches; given an 
implementation of the renormalization group they can be derived from the eigen- 
values of the linearization of the renormalization-group flow around the fixed- 
point S* in Fig. [T2] Figure [15]' a) shows our numerical estimates for several 
critical exponents of the RFIM model for Barkhausen noise [35,36], together 
with our 6 — e expansions results [34, 37]. Of course, the challenge is not to get 
analytical work to agree with numerics: it is to get theory to agree with experi- 
ment. Figure fTBT b) compares recent Barkhausen experimental measurements of 
the critical exponents to values from the three theoretical models. 



8 Notice that we can write dj and a in terms of the other exponents. These are exponent relations; 
all of the exponents can typically be written in terms of two or three basic ones. We shall derive some 
exponent relations in section [5~2l 
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5.2. Scaling functions 

Critical exponents are not everything, however. Many other scaling predictions 
are easily extracted from numerical simulations, even if they are inconvenient to 
calculate analytically. (Universality should extend even to those properties that 
we have not been able to write formulae for.) In particular, there are an abundance 
of functions of two or more variables that one can measure, which are predicted 
to take universal scaling forms. 

5.2.1. Average pulse shape 

For example, let us consider the time history V(t) of the avalanches (Fig. 0). 
Each avalanche has large fluctuations, but one can average over many avalanches 
to get a typical shape. Averaging V(t) together for large and small avalanches 
would seem silly, since only the large avalanches will last long enough to con- 
tribute at late times. The experimentalists originally addressed this issue [40] 
by rescaling the voltages of each avalanche by the peak voltage and the time by 
the duration, and then averaging these rescaled curves (dark circles in Fig.[T6b). 
But there are so many more small avalanches than large ones that it seems more 
sensible to study them separately. 

Consider the average voltage V(T, t) over all avalanches of duration T. Uni- 
versality suggests that this average should be the same for all experiments and 
(correct) theories, apart from an overall shift in time and voltage scales: 

V cxp (T,t) = AVth(T/B, t/B). (5.5) 

Comparing systems with a shifted time scale becomes simpler if we change vari- 
ables; let v(T, t/T) = V(T, t). Now, as we did for the avalanche size distribution 
in section lBTTl let us compare a system with itself, after coarse-graining by a small 
factor B = 1/(1- e): 

v{T, t/T) = Av{T/B, t/T) = (1 + a)v ((1 - e)T, t/T) . (5.6) 

Again, making e small we find av — Tdv/dT, with solution v(T,t/T) — 
vo(T)T a . Here the integration constant vq will now depend on t/T, so we arrive 
at the scaling form 

V(T,t) = T a V(t/T) (5.7) 

where the scaling function V = vq (up to an overall constant factor) is a universal 
prediction of the theory. 

Can we write the exponent a in terms of the exponents we already know? 
Since the size of an avalanche is defined as the integral of V(t), we can use 
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Fig. 16. (a) Theoretical average pulse shape scaling functions for our nucleated model and the 
front propagation model [38]. The overall height is non-universal; the two curves are otherwise 
extremely similar. The front propagation model has l/auz = 1.72 ± 0.03 in this collapse; our 
nucleation model has l/auz = 1.75 ± 0.03 (there is no reason to believe these two should agree). 
The inset shows the two curves rescaled to the same height (the overall height is a non-universal 
feature): they are quantitatively different, but far more similar to one another than either is to the 
experimental curves in part (b). (b) Comparison of experimental average pulse shapes for fixed 
pulse duration, as measured by three different groups [21,38—40]. Notice that both theory curves 
are much more symmetric than those of the experiments. Notice also that the three experiments 
do not agree. This was a serious challenge to our ideas about the universality of the dynamics of 
crackling noise [1]. (c) Pulse shape asymmetry experiment [19]. Careful experiments show a 
weak but systematic duration dependence in the collapse of the average Barkhausen pulse shape. The 
longer pulses (larger avalanches) are systematically more symmetric (approaching the theoretical 
prediction), (d) Pulse shape asymmetry theory [19]. Incorporating the time-retardation effects 
of eddy currents into the theoretical model produces a similar systematic effect. The non-universal 
effects of eddy currents are in principle irrelevant for extremely large avalanches. 



28 



J. P. Sethna 



the scaling relation (eq I5.71 i to write an expression for the average size of an 
avalanche of duration T, 



We also know that avalanches have fractal dimension df = 1 jov, so S ~ L x l av , 
and that the duration of an avalanche of size L goes as T ~ i z . Hence <S(T) ~ 
y l/cfa ^ T a+1 , and a = 1/ avz — 1. This is an example of an exponent relation. 

Can we use the experimental pulse shape to figure out which theory is cor- 
rect? Fig. fToT a) shows a scaling collapse of the pulse shapes for our RFIM and 
the front propagation model. The scaling collapse tests the scaling form eq 15.71 
by attempting to plot the scaling function V(t/T) = T~ a V(t, T) for multiple 
durations T on the same plot. The two theoretical pulse shapes look remarkably 
similar to one another, and almost perfectly time-reversal symmetric. The mean- 
field model also is time-reversal symmetric]^ Thus none of the theories describe 
the strongly skewed experimental data (Fig. [Toll. Indeed, the experiments did not 
even agree with one another, calling into question whether universality holds for 
the dynamical properties. 

This was recognized as a serious challenge to our whole theoretical picture [1, 
38,43]. An elegant, convincing physical explanation was developed by Colaiori 
et al. [44], who attribute the asymmetry to eddy currents, whose slow decay lead 
to a time-dependent damping of the domain wall mobility. Incorporating these 
eddy current effects into the model leads to a clear correspondence between their 
eddy-current theory and experiment, see Fig. fToT c.d). In those figures, notice that 
the scaling 'collapses' are imperfect, becoming more symmetric for avalanches 
of longer durations. These eddy-current effects are theoretically irrelevant per- 
turbations; under coarse-graining, they disappear. So, the original models are 
in principle correct^ but only for avalanches far larger and longer-lasting than 
those actually seen in the experiments. 

5.2.2. Avalanche size distribution 

Finally, let us conclude by analyzing the avalanche size distribution in our RFIM 
for hysteresis. The avalanche size distribution not only connects directly to many 
of the experiments (Fig. [8}, it also involves a non-obvious example of a scaling 

'The mean-field model apparently has a scaling function which is a perfect inverted parabola [41]. 
The (otherwise similar) rigid-domain-wall model, interestingly, has a different average pulse shape, 
that of one lobe of a sinusoid [42, eq. 3.102]. 

10 One might wonder how the original models do so well for the avalanche size distributions and 
other properties, when they have such problems with the pulse shape. In models which obey no 
passing [45, 46], the domains flipped in an avalanche are independent of the dynamics; eddy currents 
in these models won't change the shapes and sizes of the avalanches. 




(5.8) 
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Fig. 17. Avalanche size distribution. The distribution of avalanche sizes in our model for hysteresis. 
Notice the logarithmic scalesQ3 &) Although only at R c ~ 2.16 do we get a pure power law (dashed 
line, D(S) oc S~ T ), we have large avalanches with hundreds of domains even a factor of two away 
from the critical point, and six orders of magnitude of scaling at 5% above R c . (ii) The curves have 
the wrong slope except very close to the critical point. Be warned that a power law over two decades 
(although often publishable [47]) may not yield a reliable exponent, (iii) The scaling curves (thin 
lines) work well even far from R c . Inset: We plot D(S) / S~ r versus S' T (R — R c )/R to extract 
the universal scaling curve T>(X) (eqn l5.15l . (We use the exponent relation f = r + <t/35 without 
deriving it; see footnote|6]on page !24l ) Varying the critical exponents and R c to get a good collapse 
allows us to measure the exponents far from R c , where power-law fits are still unreliable. 
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variable). Fig. [TTI shows the distribution of avalanche sizes D[ nt (S, R) for dif- 
ferent disorders R > R c , integrated over the entire hysteresis loop (Pig. |6j. We 
observe a power-law distribution near R = R c ~ 2.16, that is cut off at smaller 
and smaller avalanches as we move away from R c 

Let us derive the scaling form for Di nt (S,R). By using scale invariance, 
we will be able to write this function of two variables as a power of one of the 
variables times a universal, one-variable function of a combined scaling variable. 
From our treatment at R c (eqns !5. lb we know that 

S' = S/(l + ce), 

' (5 9) 

D' = D{1 + ae) . 1 ' 

A system at R = R c + r after coarse-graining will have all of its avalanches 
reduced in size, and hence will appear similar to a system further from the critical 
disorder (where the cutoff in the avalanche size distribution is smaller, Fig.|7]l, say 

at R = R c + Er = R c + (1 + ee)r. Hence 

D(S', R c + Er) = D'(S', R c + r) = AD(CS', R c + r), 
D{S',R c + {l + ee)r) = (l + ae)D((l + ce)S',R c + r). °' ' 

To facilitate deriving the scaling form for multiparameter functions, it is help- 
ful to change coordinates so that all but one variable remains unchanged un- 
der coarse-graining (the scaling variables). In the average pulse shape of sec- 
tion 15.2.11 the time t and the duration T change in the same way under coarse- 
graining, so the ratio was a scaling variable. For the avalanche size distribution, 
consider the combination X = S e ' c r. After coarse-graining S' = S/C and 
shifting to the higher disorder r' = Er this combination is unchanged: 



X' = S' e/c r' = (S/C) e/c {Er) = (S/ (1 + ce)) e/c ((1 + ee)r) 

1 + ee 
Xl + ce) e ' 



= S e/C r ( 77-7— rrTT ) = S e/C r + 0(e 2 ) = X + 0( e 2 ). 



(5.11) 



Let D(S, S e / C R) = D{S, R) be the size distribution in terms of the variables S 
and X. Then D coarse-grains much like a function of one variable, since X stays 
fixed. Equation 15 . 1 Ql no w becomes 

D(S', X') = D(S', X) = (l + ae) D ((1 + ce) S' , X) , (5.12) 

so 

BD 



13 This cutoff is the one described in section \5A\ scaling as (R — R c ) " . 
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and hence 

D(S,X) = S~ a/C = S- f V(X) (5.14) 

with the universal scaling function T>(X). Rewriting things in terms of the orig- 
inal variables and the traditional Greek names for the scaling exponents (c = 
l/crv, a = f/crv, and e = 1/v), we find the scaling form for the avalanche size 
distribution: 

D(S,R) «.S~ f V(S"(R- Rc)). (5.15) 

We can use a scaling collapse of the experimental or numerical data to extract 
this universal function, by plotting D/S~ T against X = S a (R — R c ); the inset 
of Fig. [T71 shows this scaling collapse. 

In broad terms, most properties that involve large scales of length and time at 
a critical point will have universal scaling forms; any N variable function will be 
writable in terms of a power law times a universal function of N — 1 variables, 
F(x, y, z) ~ z~ a T{xl ' z 13 , y/z 1 ). The deep significance of the renormalization- 
group predictions are only feebly illustrated by the power-laws most commonly 
studied. The universal scaling functions, and other morphological self-similar 
features are the best ways to measure the critical exponents, the sharpest tests for 
the correctness of theoretical models, and provide the richest and most complete 
description of the complex behavior observed in these systems. 

References 

[1] J. P. Sethna, K. A. Dahmen, and C. R. Myers. Crackling noise. Nature, 410:242, 2001. 

[2] James P. Sethna, Karin A. Dahmen, and Olga Perkovic. Random-field models of hysteresis, 
http: //arxiv. org/abs/cond-mat/0406320 In The Science of Hysteresis, Vol. II, 
pages 107-179. Academic Press, 2006. 

[3] James P. Sethna. Statistical Mechanics: Entropy, Order Parameters, and Complex- 
ity, http : / '/www .physics . Cornell . edu/ sethna /StatMech/ Oxford University 
Press, Oxford, 2006. 

[4] M. C. Kuntz, P. Houle, and J. P. Sethna. Crackling noise, simscience . org/ 1998. 

[5] B. Gutenberg and C. F. Richter. Seismicity of the Earth and Associated Phenomena. Princeton 

University, Princeton, NJ, 1954. 
[6] P. A. Houle and J. P. Sethna. Acoustic emission from crumpling paper. Physical Review E, 

54:278, 1996. 

[7] L. I. Salminen, A. I. Tolvanen, and M. J. Alava. Acoustic emission from paper fracture. Phys. 
Rev. Lett., 89:185503, 2002. 

[8] M. Carmen Miguel, A. Vespignani, S. Zapperi, J. Weiss, and J.-R. Grasso. Intermittent dislo- 
cation flow in viscoplastic deformation. Nature, 410:667, 2001. 

[9] F. Kun, Gy. B. Lenkey, N. Takacs, and D. L. Beke. Structure of magnetic noise in dynamic 
fracture. Phys. Rev. Lett., 93:227204, 2004. 
[10] E.M. Kramer and A.E. Lobkovsky. Universal power law in the noise from a crumpled elastic 
sheet. Phys. Rev. E, 53:1465-1469, 1996. 



32 J. P. Sethna 

[11] L. C. Krysac and J. D. Maynard. Evidence for the role of propagating stress waves during 
fracture. Phys. Rev. Lett., 81:4428-31, 1998. 

[12] Knut J0rgen Mal0y, Stephane Santucci, Jean Schmittbuhl, and Renaud Toussaint. Local waiting 
time fluctuations along a randomly pinned crack front. Phys, Rev. Lett., 96:045501, 2006. 

[13] S. Tewari, D. Schiemann, D.J. Durian, CM. Knobler, S.A. Langer, and A.J. Liu. Statistics of 
shear-induced rearrangements in a two-dimensional model foam. Phys. Rev. E, 60:4385^-396, 
1999. 

[14] N. Vandewalle, J.F. Lentz, S. Dorbolo, and F. Brisbois. Avalanches of popping bubbles in 

collapsing foams. Phys. Rev. Lett, 86:179, 2001. 
[15] M. Cieplak and M.O. Robbins. Dynamical transition in quasistatic fluid invasion in porous 

media. Phys. Rev. Lett, 60:2042-2045, 1988. 
[16] O. Narayan and D.S. Fisher. Threshold critical dynamics of driven interfaces in random media. 

Phys. Rev. B, 48:7030-42, 1993. 
[17] Michael D. Uchic, Dennis M. Dimiduk, Jeffrey N. Florando, and William D. Nix. Sample 

dimensions influence strength and crystal plasticity. Science, 305:986-9, 2004. 

[18] Dennis M. Dimiduk, Chris Woodward, Richard LeSar, and Michael D. Uchic. Scale-free inter- 
mittent flow in crystal plasticity. Science, 312:1188-90, 2006. 

[19] Stefano Zapped, Claudio Castellano, Francesca Colaiori, and Gianfranco Durin. Signature of 
effective mass in crackling-noise asymmetry. Nature Physics, 1:46-9, 2005. 

[20] James P. Sethna, Karin Dahmen, Sivan Kartha, James A. Krumhansl, Bruce W. Roberts, and 
Joel D. Shore. Hysteresis and hierarchies: Dynamics of disorder-driven first-order phase trans- 
formations. Physical Review Letters, 70:3347-50, 1993. 

[21] G. Durin and S. Zapped. Scaling exponents for barkhausen avalanches in polycrystalline and 

amorphous ferromagnets. Phys. Rev. Lett., 84:4705-4708, 2000. 
[22] H. Ji and M. O. Robbins. Percolative, self-affine, and faceted domain growth in random 3- 

dimensional magnets. Phys, Rev. B, 46:14519, 1992. 
[23] O. Narayan. Self-similar barkhausen noise in magnetic domain wall motion. Phys. Rev. Lett., 

77:3855-3857, 1996. 

[24] B. Koiller and M. O. Robbins. Morphology transitions in three-dimensional domain growth 
with gaussian random fields. Phys, Rev. B, 62:5771-5778, 2000. and references therein. 

[25] B. Alessandro, C. Beatrice, G. Bertotti, and A. Montorsi. Domain-wall dynamics and 
barkhausen effect in metallic ferromagnetic materials. 2. experiments. J. Appl. Phys., 68:2908- 
2915, 1990. 

[26] B. Alessandro, C. Beatrice, G. Bertotti, and A. Montorsi. Domain-wall dynamics and 
barkhausen effect in metallic ferromagnetic materials. 1. theory. /. Appl. Phys., 68:2901-2907, 
1990. 

[27] R Cizeau, S. Zapperi, G. Durin, and H.E. Stanley. Dynamics of a ferromagnetic domain wall 

and the barkhausen effect. Phys. Rev. Lett., 79:4669, 1997. 
[28] S. Zapperi, P. Cizeau, G. Durin, and E. H. Stanley. Dynamics of a ferromagnetic domain wall: 

Avalanches, depinning transition, and the barkhausen effect. Phys. Rev. B, 58:6353-6366, 1998. 
[29] Gianfranco Durin. The barkhausen effect, [http : / / www . ien . it /~durin/barkh . html| 

2001. 

[30] E. A. Guggenheim. The principle of corresponding states. Journal of Chemical Physics, 
13:253-61, 1945. 

[31] R Heller and G. B. Benedek. Nuclear magnetic resonance in MnF2 near the critical point. 
Physical Review Letters, 8:428-32, 1962. 



Crackling Noise and Avalanches 



33 



[32] J. Zinn- Justin. Quantum field theory and critical phenomena (3rd edition). Oxford University 
Press, 1996. 

[33] M. E. J. Newman. Power laws, Pareto distributions and Zipf's law. Contemporary Physics, 

46:323-51, 2005. [arXiv.org/abs/cond-mat/04 12 004 /| 
[34] K. Dahmen and J. P. Sethna. Hysteresis, avalanches, and disorder induced critical scaling: A 

renormalization group approach. Phys. Rev. B, 53:14872, 1996. 
[35] O. Perkovic, K. Dahmen, and J. P. Sethna. Avalanches, barkhausen noise, and plain old criti- 

cality. Phys. Rev. Lett., 75:4528, 1995. 
[36] Olga Perkovic, Karin A. Dahmen, and James P. Sethna. Disorder-induced critical phenomena 

in hysteresis: Numerical scaling in three and higher dimensions. Phys. Rev. B, 59:6106-19, 

1999. 

[37] K. Dahmen and J. P. Sethna. Hysteresis loop critical exponents in 6-epsilon dimensions. Phys. 

Rev. Lett., 71:3222-5, 1993. 
[38] A. Mehta, A. C. Mills, K. A. Dahmen, and J. P. Sethna. Universal pulse shape scaling function 

and exponents: Critical test for avalanche models applied to barkhausen noise. Phys, Rev. E, 

65:046139/1-6, 2002. 

[39] L. Dante, G. Durin, A. Magni, and S. Zapped. Low-field hysteresis in disordered ferromagnets. 

Phys. Rev. B, 65:144441, 2002. 
[40] D. Spasojevic, S. Bukvic, S. Milosevic, and H. E. Stanley. Barkhausen noise: Elementary 

signals, power laws, and scaling relations. Phys. Rev. E, 54:2531-2546, 1996. 
[41] M. Kuntz. Barkhausen noise: simulations, experiments, power spectra, and two dimensional 

scaling. PhD thesis, Cornell University, 1999. 
[42] G. Durin and S. Zapperi. The barkhausen effect. In G. Bertotti and I. Mayergoyz, editors, 

Science of Hysteresis, Vol. II. Elsevier, London, 2004. 
[43] M. Kuntz and J. P. Sethna. Noise in disordered systems: The power spectrum and dynamics 

exponents in avalanche models. Phys. Rev. B, 62:11699, 2000. 
[44] F. Colaiori, M. J. Alava, G. Durin, A. Magni, and S. Zapperi. Phase transitions in a disordered 

system in and out of equilibrium. Phys, Rev. Lett., 92:257203, 2004. 
[45] Alan A. Middleton. Asymptotic uniqueness of the sliding state for charge-density waves. Phys- 
ical Review Letters, 68:670, 1992. 
[46] A. A. Middleton and D.S. Fisher. Critical behavior of charge-density waves below threshold: 

numerical and scaling analysis. Phys, Rev. B, 47:3530-52, 1993. 
[47] O. Malcai, D. A. Lidar, 0. Biham, and D. Avnir. Scaling range and cutoffs in empirical fractals. 

Physical Review E, 56:2817-28, 1997. 



